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ABSTRACT 

We study the constraints on reionization from five years of WMAP data, parametrizing the cvohition 
of the average fraction of ionized hydrogen with principal components that provide a complete basis 
for describing the effects of reionization on large-scale _E-mode polarization. Using Markov Chain 
Monte Carlo methods, we find that the resulting model-independent estimate of the total optical 
depth is nearly twice as well determined as the estimate from 3-year WMAP data, in agreement with 
simpler analyses that assume instantaneous reionization. The mean value of the optical depth from 
principal components is slightly larger than the instantaneous value; we find r = 0.097 ± 0.017 using 
only large-scale polarization, and r = 0.101 ± 0.019 when temperature data is included. Likewise, 
scale invariant n^ = 1 spectra are no longer strongly disfavored by WMAP alone. Higher moments of 
the ionization history show less improvement in the 5-year data than the optical depth. By plotting 
the distribution of polarization power for models from the MCMC analysis, we show that extracting 
most of the remaining information about the shape of the reionization history from the CMB requires 
better measurements of E-xnode polarization on scales of ^ '-^ 10 — 20. Conversely, the quadrupole 
and octopole polarization power is already predicted to better than cosmic variance given any allowed 
ionization history at z < 30 so that more precise measurements will test the ACDM paradigm. 
Subject headings: cosmic microwave background — cosmology: theory — large-scale structure of 
universe 



1. INTRODUCTION 

The amplitude of fluctuations in the _B-mode compo- 
nent of cosmic microwave background (CMB) polariza- 
C^ . tion on large scales provides the current best constraint 
^^ ' on the Thomson scattering optical depth to reionization, 
Cn I T. Assu ming that the univer se was reionized instanta- 
■»^ , neously, iDunklev et al.l (|2008[ ) estimate the total optical 
(3 ■ depth to be r = 0.087 ± 0.017 using five years of WMAP 
OO ', data. Theoretical studies suggest that the process of 
f^ . reionization was too c omplex to be well desc ribed as a 
L| ' sudden transition (e.g.. lBarkana fc Loebll200ll ). Previous 
. S^ \ studies have examined how the constraint on t depends 
\^ , on the evolution of the globally-averaged ionized frac- 
Vh ' tion during reionization, Xe{z)^ for a variety of specific 
C^ \ theoretical scenarios. If the assumed form of Xe{z) is in- 
correct, the estimated value of r can be biased; this bias 
can be lessened by considering a wider variety of reion- 
ization histories at the expense of increasing the uncer- 
taintv in r (iKa plinghat et al. 2003: Holder et al. 200^; 
[CoTombo et 311120051 ). 

For the previous release of three years of WMAP data, 
using a more model-independent analysis of reioniza- 
tion based on principal components does not change the 
basic conclusions of studies that assume instantaneous 
reionization or other simple models; however, as the po- 
larization power spectrum becomes better determined, 
it is increasingly important to adopt an approach with 
sufficient freedom to approximate a variety of possible 
reion ization scenarios in order to minimize parameter bi- 
ases (jMortonson fc Hull2008al rbl). In this letter, we study 
how model-independent constraints on reionization have 
changed with the addition of two more years of data from 
WMAP and on what scales further measurements would 



have the largest impact on ionization constraints. We 
review the principal component parametrization of the 
ionization history in § [2l We then present the current 
constraints on the principal component amplitudes and 
the optical depth to reionization (§ [3]) and the range of 
polarization power spectra for general models that are 
currently allowed by the data (§ [4|). We discuss these 
results in §[5l 

2. IONIZATION PRINCIPAL COMPONENTS 

We parametrize the reionization history as a free func- 
tion of redshift by decomposing Xe{z) into its principal 
components (PCs) with respect to the Ji^-mode polariza- 
tion o f the CMB (iHu fc Holder! [2001 iMortonson fclfal 
l2008aD : 



Xe{z) 



„fid 



(z) + ^m^S'^(z), 



(1) 



where the principal components, S^{z), are the eigen- 
functions of the Fisher matrix that describes the depen- 
dence of Cf^ on Xe{z), m^ are the amplitudes of the 
principal components for a particular reionization his- 
tory, and x^'^{z) is the fiducial model at which the Fisher 
matrix is computed. The components are rank ordered 
by their Fisher-estimated variances. The lowest- variance 
eigenmode (/x = 1) is an average of the ionized frac- 
tion over the entire redshift range, weighted toward high 
z. The ^ = 2 mode measures the difference between 
the amount of ionization at high z and at low z, and 
higher modes follow this pattern with weighted averages 
of Xe{z) that oscillate with higher and higher frequency 
in redshift. The main advantage of using principal com- 
ponents as a basis for Xe{z) is that only a small number of 
the components are required to completely describe the 
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Fig. 1. — Marginalized 2D 68% and 95% CL contours for the optical depth to rcionization (r) and the amplitudes of the 5 lowest- variance 
principal components oi Xe{z) (m^, /j, = 1 — 5). Panels along the diagonals show the ID posterior probability distributions. Constraints are 
plotted for both 3-year {red dotted lines) and 5-year {blue shading, solid lines) WMAP data. In the left plot, only the low-f reionization 
peak in the B-mode polarization power spectrum is used for parameter constraints, and all parameters besides the 5 PC amplitudes and r 
are held fixed. For the constraints in the right plot, we use both temperature and polarization data and allow five additional parameters to 
vary: Cjh^, Clch^, 9a, ^a, and n^. The plot boundaries for the PC amplitudes correspond to physicality priors that exclude models that 
are unable to satisfy < Xe < 1 for any combination of the higher-variance {fi > 6) PCs. 



effects of reionization on large-scale CMB polarization, 
so we obtain a very general parametrization of the reion- 
ization history at the expense of only a few additional 
parameters. 

The principal components are defined over a limited 
range in redshift, Zmin < z < Zmax, with Xe = at 

Z > ^max and Xe = 1 at Z < Zmin- We take Zmin = 6, 

since the absence of Gunn-Peterson absorption in the 
spectra of quasars at z < 6 indicates th at the universe i s 
nearly fully ionized at lower redshifts ()Fan et al.ll2006D . 
In the MCMC analysis presented here, we always use the 
five lowest-variance principal components of Xe{z) with 
Zma,x = 30, constructed around a constant fiducial model 
of Xg'^(z) = 0.15. The amplitudes of these components 
then serve to parametrize general reionization histories 
in the a nalysis of CMB po l arizati on data. We refer the 
reader to lMortonson fc Hul (J2008a[ ) for further discussion 
of these choices and the demonstration that five compo- 
nents suffice to describe the i?-mode spectrum to better 
than cosmic variance precision. 

We impose priors on the principal component ampli- 
tudes corresponding to physical values of the ionized frac- 
ti on, < Xe < 1, acc o rding t o the conservative approach 
of iMortonson fc Hul (|2008af) . All excluded models are 
unphysical, but the models we retain are not necessarily 
strictly physical. Finally, we neglect helium reionization, 
which is a small correction at the current level of preci- 
sion but will be more impor tant for future analyses (e.g., 
[Colombo fc PierDaolil[200l . 

3. OPTICAL DEPTH CONSTRAINTS 

We examine the implications of the WMAP 5-year 
data for general models of reionization parametrized by 
principal components using a Markov Chain Monte Carlo 
analysis that mirrors our previous study of the 3-ycar 



data in IMortonson fc Hul ()2008af ). We consider con- 
straints from either large-scale polarization alone, with 
parameters that do not directly affect reionization fixed 
to values that fit the temperature data ("EE"), or from 
the full set of temperature and polarization data, varying 
the parameters of the "vanilla" ACDM model (baryon 
density i^bh^, cold dark matter density i^ch^, acous- 
tic scale 9a, scalar amplitude Ag, and scalar spectral 
tilt Us) in addition to the reionization PC amplitudes 
("TT-hTE-hEE"). In both cases, the total optical depth 
to reionization, r, is a derived parameter. 

The MCMC constraints on principal component am- 
plitudes and the derived optical depth for both of these 
cases are plotted in Fig. |T1 along with the previous 
constr aints from 3-year WMAP data (jMortonson fc Hul 
l2008al) . While there are some improvements in all 5 of 
the individual components when considering EE alone, 
these changes are not as large as the improvement in the 
optical depth constraint when all of the data are consid- 
ered. Adding both temperature data and extra parame- 
ters in going from EE only to TT+TE-f EE has the net 
effect of slightly strengthening constraints on the higher 
ranked PC amplitudes, although there is very little ef- 
fect on r. The additional constraining power for both 
3-year and 5-year data comes mainly from the measured 
temperature power spectrum at £ ~ 10 — 100, which ex- 
cludes models with additional Dopplcr effect contribu- 
tions due to narrow fea tures in the ionization history 
(jMortonson fc Hull2008a[) . 

Modeling reionizat ion as an instan t aneou s transition at 
some redshift Zrcion, iDunklev et al.l ()2008[ ) estimate the 
optical depth from the 5-year WMAP data to be r = 
0.087±0.017, almost a factor of two more precise than th e 
estimate from three years of data (jSpergel et al.ll2007l ). 
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Fig. 2. — Median {blue solid curve) and 68% and 95% CL regions [blue shaded regions) of polarization power spectra for any ionization 
history at z < 30 allowed by the 5-year WMAP data, assuming the standard ACDM paradigm. As in Fig. [T] only polarization data are 
used in the left panel and both temperature and polarization data are included in the right panel. Red dashed curves: cosmic variance 
(68% and 95% CL) around the median model. 



For ionization histories parametrized by PCs, we find 
that the constraint on optical depth is r = 0.097 ±0.017 
for the EE case and r = 0.101 ± 0.019 for TT+TE+EE. 
As with the 3-year data, the error on r is roughly 10% 
larger with the inclusion of temperature data and a larger 
set of parameters, and the error in both cases is the same 
or only slightly larger than for the instantaneous reion- 
ization analysis. 

The central value of r for more general ionization 
histories is higher than the instantaneous reionization 
value by~ 0.5 — Icr; a similar shift toward larger 
optica l depths was seen in P C analysis of the 3-year 
data (jMortonson fc Hul l2008a| ) . The maximum likeli- 
hood model, however, has r = 0.088 for EE data and 
T = 0.090 for TT+TE+EE, much closer to the instanta- 
neous reio nization maxirnum li kelihood optical depth of 
r = 0.089 (jDunklev et al.ll2008f l. The larger mean optical 
depth is at least partly due to having a large parameter 
volume of models with finite ionization fraction at high 
redshift that are still allowed by the data. With flat pri- 
ors on the principal components, this volume effect can 
boost the mean optical depth of models even though the 
mean likelihood of low optical depth models remains the 
same. Our assumption of full ionization at redshifts be- 
low Zmin = 6 for all models also limits how small the 
optical depth can be. Relative to the best-fit instan- 
tan eous model with Zreion = 11-0 ± 1.4 (jDunklev et al.l 
l2008i) . there are simply more ways to increase r than 
there are to decrease r by changing the ionization his- 
tory, given these priors and the current data. 

For the TT+TE+EE analysis, the larger mean optical 
depth is accompanied by shifts in correlated parameters, 
particularly the spectral tilt: n^ = 0.990 ± 0.024 with 
Xe{z) parametrized by P Cs, and Ug = 0.960+0. 015 for in- 
stantaneous reionization (jKomatsu et al.|[2008( l. As with 
the optical depth, however, some of this shift is a parame- 
ter volume effect. The maximum likelihood model for the 
principal component analysis has Us = 0.976. The best 
fit scale invariant model (fixing Ug = 1) is a poorer fit 



to the data by Ax|ff = -21n(/://:i„ax) ~ 1, where £,„ax 
is the maximum likelihood. (The instantaneous reion- 
ization maximum likelihood model is also at Ax^^j w 1 
relative to the best fit with principal components.) As 
measurements of CMB polarization improve with future 
data, particularly with detections in the 10 < ^ < 20 
range (see Sj4l): the constraints on parameters such as op- 
tical depth and tilt should become less sensitive to our 
assumptions about the priors. 

Unlike the total optical depth, constraints on the op- 
tical depth over more limited redshift ranges have only 
improved slightly. With three years of WMAP data, the 
95% upper limit on the optical depth from 2; > 20 (al- 
lowing for a significant ionized fraction iip to z ~ 40) 
was t{z > 20) < 0.08 (jMortonson fc Hull2008a[) . The 
limit from 5-year data is t{z > 20) < 0.07. If we instead 
choose the dividing redshift to be the best-fit value of 
the redshift of instantaneous reionization, Zrcion = 11, we 
find a similar constraint for the contribution to the op- 
tical depth from high redshift: t{z > 11) < 0.07. Com- 
pared to the 3-year data, there is also a more significant 
(but still weak) preference for nonzero optical depth from 
6 < z< 11. 

4. POWER SPECTRA OF ALLOWED MODELS 

To better understand at which scales the reionization 
peak of £'-mode polarization is best constrained by the 
current data, we plot the 68% and 95% CL limits on 
Cf^ from the Monte Carlo chains in Fig. [2] These limits 
reflect the range of ensemble-averaged power allowed by 
the 5-ycar data and the PC-parametrized reionization 
histories. Since this parametrization is complete in the 
power spectrum, the range in Fig. |2| reflects the allowed 
model power spectra for any ionization history at z < 

At ^ < 5, the variation in allowed models is smaller 
than the uncertainty due to cosmic variance. In other 
words, the data at £ ^ 5 in combination with any ion- 
ization history and the power law initial power spectrum 



make a prediction for the ensemble-averaged power at 
lower i that is sharper than can be measured. Conversely, 
measurements that violate this prediction at a statisti- 
cally significant level require modifications to the ACDM 
paradigm itself, much like low measurements of the tem- 
perature quadrupole. It is interesting that the maximum 
likelihood _E-mode polarization quadrupole reported by 
iNolta et all (l200l . 6Cf -^7277 w 0.15 /zK^, is in excess of 
the 95% cosmic variance region shown in Fig. [21 

The uncertainty in the model space is largest at in- 
termediate scales of £ ~ 10 — 20, where the large-scale 
polarization power is expected to be smallest. There is 
substantial room for improved measurements of the spec- 
trum on these scales before reaching the cosmic vari- 
ance limit. Tighter constraints on the iJ-mode power 
at 10 < ^ < 20 would better determine the amplitude 
of principal components beyond the first; such measure- 
ments are necessary to be able to discriminate among dif- 
ferent reionization histories with the same total optical 
depth. Physically, these measurements would better con- 
strain the ionization history at high redshifts (z > 15). 

On small scales (£ > 30 — 40) the limits on power spec- 
tra of reionization models from the chains again become 
tighter than cosmic variance since the theoretical ampli- 
tude of the recombination peak is well determined due to 
constraints on parameters from the temperature data. 

Comparison of the two panels in Fig. [2] shows that the 
main effect of including temperature data in constraints 
on principal components is to eliminate models with large 
power at 10 < £ < 30. As mentioned in the previous sec- 
tion, these models are excluded by the data due to their 
increased temperature fluctuations at £ ~ 10—100. Even 
with TT-I-TE-I-EE data, the range of power in models al- 
lowed by current data is a few times larger than cosmic 
variance. 



5. DISCUSSION 

The 5-year WMAP polarization data significantly im- 
prove the estimate of the total optical depth, reducing 
the error from Cr ~ 0.03 to CTt « 0.017. This improve- 
ment is seen in both a model-independent analysis us- 
ing principal components of the ionization history and in 
an analysis that assumes instantaneous reionization, al- 
though there is a small shift in the central value with the 
model-independent method preferring a slightly higher 
mean around t = 0.1. 

As with the 3-year data, the £^-mode reionization peak 
is currently best measured on the largest scales, £ ~ 5. 
Determining details of the ionization history beyond the 
optical depth requires information about the full shape of 
the reionization peak, which can be obtained by supple- 
menting the current observations with better measure- 
ments of the E-mode power on scales of 5 < ^ < 30. In 
particular, improved knowledge of the power on scales 
between the main reionization peak and the recombina- 
tion peak dX I ^ 10 — 20 would be the most useful for 
distinguishing models of reionization with different ion- 
ization histories but the same optical depth. 

Conversely, the data along with any allowed ioniza- 
tion history at z < 30 in the standard ACDM con- 
text already predict the ensemble-averaged polarization 
quadrupole and octopole powers to better than cosmic 
variance. More precise measurements in this regime can 
test the standard model itself. 
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